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We report on previously unobserved inter-system degeneracies in two- component equal-mass Fermi 
gases with interspecies zero-range interactions under isotropic harmonic confinement. Over the 
past 10 years, two-component Fermi gases consisting of n± spin-up and 712 spin-down atoms with 
interspecies zero-range interactions have become a paradigm for modeling condensed matter systems, 
nuclear matter and neutron matter. We show that the eigen energies of the (m + 1, ni — 1) system 
are degenerate with the eigen energies of the (711,712) system for any s-wave scattering length a s , 
including infinitely large, positive and negative a s . The existence of the inter-system degeneracies 
is demonstrated explicitly for few-body systems with m +712 = 4, 5 and 6. The degeneracies and 
associated symmetries are explained within a group theoretical framework. 
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Symmetry is one of the most fundamental concepts in 
physics, underlying our understanding of elementary par- 
ticle physics, relativity and quantum mechanics, to name 
a few [1]. In quantum mechanics, symmetries manifest 
themselves in degeneracies of energy eigen values. If a 
Hamiltonian is invariant under rotations, for example, 
the eigen energies are (2L + l)-fold degenerate, where L 
denotes the orbital angular momentum quantum num- 
ber [![. Similarly, the fact that the energy spectrum 
of the non-relativistic hydrogen atom depends only on 
the principal quantum number is intimately related to 
conserved quantities associated with the orbital angular 
momentum and Runge-Lenz vectors [2[. 

Dilute atomic two-component Fermi gases with short- 
range interspecies s-wave interactions can nowadays be 
realized routinely in many cold atom laboratories [3|- In 
these experiments, the atoms occupy two different hyper- 
fine states that are interpreted as spin- 1/2 pseudo states. 
Ultracold atomic Fermi gases have emerged as model sys- 
tems with which to study condensed matter phenomena 
such as the BCS-BEC crossover and nuclear physics phe- 
nomena such as the equation of state of superfluid neu- 
tron matter jj-[6|. A multitude of results have been ob- 
tained for two-component equal-mass Fermi gases with s- 
wave zero-range (ZR) interactions. A notable milestone 
is the derivation of various universal relations by Tan (Ti- 
ll)], which are centered around the "contact" and now 
form the basis for novel spectroscopic techniques [Io|, Ell • 
Another notable milestone is the identification of a hid- 
den SO (2,1) symmetry of the two-component Fermi gas 
with ZR interactions at unitary in an isotropic harmonic 
trap by Werner and Castin [l2[, which manifests itself in 
ladders of uniformly spaced excitation frequencies. 

Our work identifies another symmetry that manifests 
itself in the existence of degenerate eigen energies of two- 
component equal-mass Fermi gases with the same num- 
ber of particles but different numbers of spin-up and spin- 
down atoms, i.e., of (771,772) and (n'^n^) systems with 



emerge in the ZR limit for any value of the interspecies 
5-wave scattering length and are broken for finite-range 
interactions or unequal-mass systems. 

Our starting point is the non-relativistic Hamiltonian 
H of the two-component Fermi gas with n 1 spin-up atoms 
and 772 spin-down atoms (77 = 77 1 + 772), 



H = H + V in 
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and Vi n t describes the interactions between the spin-up 
and spin-down atoms, 
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In Eq. (pQ), 777 denotes the atom mass, uo the angular 
trapping frequency, and fj the position vector of the jth 
atom measured with respect to the trap center. Follow- 
ing the literature [J], the spin- up and spin-down compo- 
nents by themselves are assumed to be non-interacting 
(NI). We model the intercomponent atom-atom interac- 
tions by a short-range Gaussian potential V g [13] with 
depth V and range r , V g (r jk ) = -V exp[-r 2 /(2r^)], 
where rjk = \rj — rk\- For a fixed 7*0, we adjust the 
depth Vo such that V g reproduces the desired free-space 
zero-energy atom-atom s- wave scattering length a s . We 
restrict ourselves to two-body potentials that support no 
two-body 5-wave bound state in free-space for negative 
a s and one two-body s- wave bound state in free-space for 
positive a s . In the ro —> limit, our interaction model 
provides a realization of the ZR ^-function interaction. 
In practice, we determine the eigen energies of H for a 
sequence of ro values and then extrapolate the eigen ener- 
gies to the ro —> limit. Throughout, we consider ranges 
ro that are much smaller than the harmonic oscillator 
length a no , where a no = \/h/(muj). 
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FIG. 1: (Color online) Illustration of the n — 4 systems. 
Panels (a) and (b) show the ground state configurations of 
the NI trapped (3, 1) and (2, 2) systems. The horizontal solid 
lines indicate the single particle harmonic oscillator orbitals 
with energy (2n + I + 3/2)^j; the (n, Z) = (0,0) and (0, 1) 
orbitals are respectively one-fold and three-fold degenerate. 
Solid lines in panels (c) and (d) illustrate the spin-up — spin- 
down interactions of the (3, 1) and (2, 2) systems. 



We first consider the Hamiltonian Hq, which describes 
n NI particles under isotropic harmonic confinement. As 
an example, Figs. QJa) and DJb) illustrate the ground 
state configurations of the (77-1,77-2) = (3,1) and (2,2) 
systems. The lowest single particle orbital with energy 
3hu/2 can be occupied by a spin- up atom and a spin- 
down atom. To obey the Pauli exclusion principle, the 
other spin-up and spin-down atoms need to occupy one of 
the three excited state orbitals with energy 5hu/2. This 
simple picture yields a ground state energy of 8hu for 
both the (3, 1) and (2, 2) systems. It can be readily shown 
that the ground state of the (3, 1) system has L u = 1 + 
symmetry, where L denotes the orbital angular momen- 
tum quantum number and II the parity; this ground state 
is three-fold degenerate due to the rotational invariance 
of the Hamiltonian. The ground state of the (2, 2) system 
is nine- fold degenerate [14[. Just as the NI ground state 
manifolds of the (3, 1) and (2, 2) systems contain degen- 
erate energies corresponding to the same L u symmetry, 
so do the NI excited state manifolds. Moreover, anal- 
ogous degeneracies are readily identified for NI systems 
with larger n. 

In this paper, we are interested in the "inter- system 
degeneracies", i.e., in the fact that the (77-1 -f 1,77-2 — 1) 
and (77-1,77-2) systems support degenerate energies corre- 
sponding to the same L u symmetry. Specifically, we 
analyze what happens to the inter-system degeneracies 
when the interactions are turned on. For example, since 
the (3,1) system contains three spin- up — spin-down pairs 
while the (2, 2) system contains four [see Figs. [TJc) and 
QJd)], it seems natural to expect that the interactions 
break the inter-system degeneracies discussed above for 
the NI 77- = 4 systems. As we will show, however, this is 
not the case if ro is taken to zero: For ZR interactions, 
the eigen energies of the (3, 1) system form, within our 
numerical accuracy, a subset of the eigen energies of the 




FIG. 2: (Color online) Relative extrapolated ZR energies for 
the n = 4 systems described by H as a function of aho/a s . 
Solid and dashed lines show the extrapolated ZR energies of 
the energetically lowest-lying and second lowest-lying states of 
the (3, 1) system with 1 + symmetry while squares and circles 
show those of the (2, 2) system. The energies of the (3, 1) 
and (2, 2) systems are indistinguishable on the scale shown. 
The inset shows the fractional difference Ae, Ae 
E™{)/E™{, for the energetically lowest-lying state. 



(E r 2% 



(2, 2) system. Analogous results are found for systems 
with 77- = 5 and 6. 

To determine the eigen energies of the Hamiltonian H 
for finite depth Vb of the Gaussian model potential V^, 
we resort to the stochastic variational approach (l3l flif . 
We separate the center of mass motion and expand the 
relative eigen functions in terms of a basis set with good 
orbital an gul ar momentum quantum number L and par- 
ity n [TBI- 17]. The proper permutation symmetry under 
the exchange of identical fermions is imposed by applying 
an appropriately chosen anti-symmetrization operator to 
the basis functions. Our implementation [l8| allows for 
the treatment of states with all L u symmetries. The 
stochastic variational approach results in variational up- 
per bounds to the exact eigen energies [15[. 

As an example, Fig. [2] shows the extrapolated ZR en- 
ergies for the (3, 1) and (2, 2) systems with 1 + symmetry 
as a function of the inverse s- wave scattering length a~ x . 
In this representation, the weakly-attractive BCS regime 
(a s < and \a s \/a\ lo <C 1) is realized on the left of the 
graph and the repulsive BEC regime (a s > 0) on the 
right of the graph. Lines show the relative eigen energies 
of the (3, 1) system and symbols those of the (2, 2) system 
for ro = 0. The extrapolated ZR energies are estimated 
to have a combined basis set and ZR extrapolation error 
smaller than 0.001 % for the energetically lowest-lying 
state and smaller than 0.01 % for the energetically sec- 
ond lowest-lying state. On the scale shown, the eigen 
energies of the (3, 1) and (2, 2) systems are, somewhat 
surprisingly, indistinguishable for both the lowest and 
second lowest states for all scattering lengths considered. 
The inset shows that the fractional difference is smaller 
than 2 x 10 -4 % for the energetically lowest lying 1 + state 
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FIG. 3: (Color online) Analysis of the extrapolated ZR en- 
ergies of the n = 4 systems described by H at unitarity. 
Symbols show the fractional difference As, As = (s^'J ' n — 



u,L,U\ I u,L,Il 



)/ s 3,'i ' 5 between s^'2 ' and s^ ' as a function of . 
As is of the order of or smaller than the fractional numerical 
uncertainty of the extrapolated ZR energies. 



and the scattering lengths considered. Thus, within our 
numerical accuracy, the energy curves agree throughout 
the crossover. 

To see if the (3,1) and (2, 2) energies are also degener- 
ate for other symmetries and for higher- lying excitations, 
we focus on the infinite scattering length regime. We ana- 
lyze the extrapolated ZR energies of the (3,1) and (2, 2) 
systems at unitarity for all states with relative energy 
E rel equal to or smaller than 21fkj/2, which were deter- 
mined in Ref. [HI with an accuracy of 0.1 % or better. 
In this energy window, there exist 164 and 286 eigen en- 
ergies for the (3,1) and (2,2) systems, respectively [191 ]. 
As pointed out by Werner and Castin 12], the existence 
of a hidden SO (2,1) symmetry leads to ladders of ener- 
gies spaced by 2fkj, i.e., the relative eigen energies at 
unitarity can be written as E rel = (s u,L,u + 2q + l)huj, 



where q = 0,1, 



The separation constants s u ' L,u 



arise when solving the (m, n2)-fermion problem within 
the hyperspherical framework. We find that the rela- 
tive eigen energies with E rel < 21fio;/2, corresponding 
to (3, 1) and (2, 2) states that are affected by the inter- 
actions, are characterized by 89 and 170 s u,L,u values, 
respectively [20|. Quite surprisingly, every s UjLjU value 
of the (3, 1) system, within the numerical accuracy [181 ]. 
appears in the sequence of s u,L,u values of the (2, 2) sys- 
tem. Figure [3] shows that the fractional difference be- 
tween the s u,L,u values of the (3, 1) and (2, 2) systems is 
of the order of or smaller than the numerical accuracy 
of the extrapolated ZR energies. This suggests that the 
exact ZR energies of the (3,1) system at unitarity form 
a subset of the exact ZR energies of the (2,2) system 
at unitarity. These findings are corrobated by extensive 
perturbative calculations [21| . 

The calculations presented so far strongly suggest that 
the (3, 1) energies are degenerate with a subset of the 
(2,2) energies in the ro — >• limit for all a s . The sup- 
plemental material [2l[ shows, using the stochastic varia- 
tional and perturbative approaches, that analogous inter- 
system degeneracies exist for systems with n = 5 and 6. 
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FIG. 4: (Color online) Comparison of the eigen energies of 
the Hamiltonian H and H' with Gaussian model interactions 
as a function of ro / a^o at unitarity. Crosses and circles show 
the relative eigen energies E™^ n2 of the energetically lowest- 
lying state with 1 + symmetry described by H while diamonds 
and stars show the relative eigen energies of the energetically 
lowest-lying state with 1 + symmetry described by H r (on the 
scale shown, the diamonds and stars are indistinguishable). 
Solid lines show fits to the relative eigen energies of H. 



To interpret our observations, we construct a new Hamil- 
tonian H' 



H' = H + V^ nt , 



(4) 



that reproduces the eigen energies of the {n\ + 1, 77,2 — 1) 
and (724,712) systems described by H when ro — » 0. The 
interaction potential V{ nt includes interactions between 
all atom pairs and not just between the spin- up and spin- 
down pairs, 



(5) 



j<k 



The Hamiltonian H' treats all atom pairs on equal foot- 
ing. In particular, V{ nt is the same for the (n\ + 1, n<i — 1) 
and (rii, 77,2) systems. Intuitively, it is clear that the anti- 
symmetry of the eigen functions under the exchange of 
like atoms "turns off" the interactions between the like 
atoms when ro — » 0, thereby ensuring that the energy 
spectra of H and H' are identical when ro — » 0. This 
behavior is illustrated exemplarily in Fig. 3] for the ener- 
getically lowest-lying state of the n = 4 systems with 1 + 
symmetry interacting through the Gaussian model po- 
tential. In the ro —> limit, the eigen energies of H and 
H' agree for the (2, 2) and (3, 1) systems; moreover, as 
already pointed out above, the eigen energies of the (3, 1) 
and (2,2) systems agree. For finite ro, the eigen energies 
of the (3, 1) and (2, 2) systems described by H are char- 
acterized by different slopes while the eigen energies of 
the (3, 1) and (2, 2) systems described by H' agree within 
our numerical accuracy for all ro- 

The key motivation for introducing the Hamiltonian 
H' is that it describes all n-particle systems with ZR 
8-wave interactions, regardless of the particle statistics. 
The fact that the (n\ + 1,77-2 — 1) and (77,1,77,2) systems 
are described by the same Hamiltonian allows us to tie 
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the evidenced degeneracy of the eigen energies to the ex- 
istence of a symmetry. In particular, according to quan- 
tum mechanics Q , the existence of degenerate eigen ener- 
gies of a Hamiltonian is a manifestation of an underlying 
symmetry. Since the Hamiltonian H' is invariant under 
the permutation of any pair of atoms, the inter-system 
degeneracies are intimately related to the structure of 
the permutation group S n . Group theoretical tools are 
widely used in quantum chemistry and molecular physics 
to (anti-) symmetrize the wave functions associated with 
the electronic and nuclear degrees of freedom [22|. Here, 
they are employed to analyze the properties of the Hamil- 
tonian H\ which has been shown to reproduce the eigen 
spectrum of the original Hamiltonian H. 

The Hilbert space of the (711,712) system is spanned 
by the direct product of the Hilbert spaces of the two 
single components or, in terms of Young tableaux, [l ni ] 
[V 



23|. Here, [l ni ] = [1, 1, • • • , 1] indicates the 
fully anti-symmetric tableau of the n\ spin-up fermions. 
The direct product can be decomposed into a direct sum 
of Young tableaux that consist of at most two columns 
K >n 2 ) 0, 

[i ni ] ® [i U2 ] = [i ni+n2 ] © [2, r i+n2 " 2 ] © 

[2 2 , l n i+^-4] e ... [ 2 n2 , l ni+n2 - 2n2 ]. (6) 

If we replace rt\ and in Eq. ^ by ri\ + 1 and — 1, re- 
spectively, and then compare with the decomposition for 
the (711,712) state space, we find that the decomposition 
of the (711,712) state space contains the decomposition of 
the (n\ + 1, 712 — 1) state space, 

[i ni ] © [i U2 ] = ([i ni+1 ] © [r 2 - 1 ]) © [2 n2 , i ni - n2 ] . (7) 

This decomposition into irreducible representations 
shows explicitly that the decomposition of the (711 + 
l,ri2 — 1) system is contained in that of the (711,712) 
system. Correspondingly, the eigen energies of the 
(71 1 + 1, 712 — 1) system with ZR interactions form a subset 
of those of the (711,712) system with ZR interactions for 
all a s . Equation ([7]) shows, in agreement with our earlier 
discussion, that the (711,712) system contains additional 
eigen energies. 

In summary, we have identified and interpreted inter- 
system degeneracies of two-component Fermi gases with 
ZR interactions under spherically symmetric confine- 
ment. The fact that the eigen energies of the 71 = 4 sys- 
tem with spin projection quantum number Ms = 1 form 
a subset of the eigen energies of the 71 = 4 system with 
Ms = (and similarly for 71 > 4) has multiple implica- 
tions. From a computational point of view, the degenera- 
cies can be used to test the accuracy of various schemes 
employed to solve the 71-fermion Schrodinger equation. 
Moreover, in certain cases it may be easier to treat the 
energetically lowest lying state of a system with larger 
Ms than an excited state of a system with smaller Ms, 
allowing one to substitute an excited state calculation by 



a ground state calculation for a system of the same size 
but with different Ms- The inter-system degeneracies 
also have experimentally observable implications. Since 
the change of the energy with scattering length coincides 
for certain eigen states of the (711 — 1, 712 + 1) and (711, 712) 
systems, the corresponding eigen states, which charac- 
terize two distinctly different physical systems, have the 
same contact [7H9|. Moreover, the two distinctly dif- 
ferent systems share a common set of eigen frequencies. 
These frequencies can be measured via microwave spec- 
troscopy [26|. 

The discussed inter-system degeneracies do not only 
exist for systems with ZR interactions but also for sys- 
tems with finite-range interactions such as electronic sys- 
tems, provided the Hamiltonian under study is invariant 
under permutation of all particle pairs. This also im- 
plies that the degeneracies are not limited to harmoni- 
cally confined systems but also exist for systems in free 
space or under non-harmonic confinement, provided the 
Hamiltonian under study is invariant under permutation. 
We conclude by noting that our analysis is based on the 
assumption that the interaction potential is constructed 
from pairwise two-body interactions. The presence of 
three-body forces, which are needed to describe non- 
universal Efimov states or nuclear systems, introduces 
a new degree of freedom not considered here. 

Support by the NSF through grant PHY-0855332 and 
by the ARO is gratefully acknowledged. 
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This material summarizes details of our calculations 
and results that support the conclusions drawn in the 
main part of the paper. Sections U and El show, within 
a perturbative framework, that the inter-system degen- 
eracies exist in the weakly-interacting limits. Section ITTT1 
provides numerical examples for the inter-system degen- 
eracies for two-component gases at unitarity with n = 4 
and 5 atoms. Last, Sec. IIVI discusses the application of 
the group theoretical framework discussed in the main 
text to systems with up to n = 6 atoms. 



I. WEAKLY- ATTRACTIVE FERMI GAS 

When |a s |/aho <^ 1 and a s < 0, the two-component 
Fermi gas can be described perturbatively assuming 5- 
function interactions between the spin-up and spin-down 
atoms. For the perturbative analysis of this regime, we 
assume that the two-body potential Vtb(^fe) m Eq. (3) 
of the main text is written as (Airh 2 ] a s / 'm)S(fjk) • Within 
first-order degenerate perturbation theory, the relative 
energy E^ of the (724,712) system can be written in 
terms of the non-interacting energy E^}^\ and an energy 
shift that is proportional to a s /ah Q QUll; 



nrel ,: 



77>rel __ Ein 

rt\ ,ri2 ni,n2 



(27T) 



"1/2, 
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The coefficients c ni ^ n2 are system and state dependent. 
The non-interacting wave functions can be constructed 
within the hyperspherical coordinate approach [3] or by 
writing the non-interacting wave functions in terms of a 
product of determinants built using the single particle 
harmonic oscillator orbitals p], Q. We checked explic- 
itly for n = 4 — 6 and all possible (711,712) combinations 
that the energy shifts are the same for the interaction 
potentials Vint and V^ nt , supporting our finding that the 
Hamiltonians H and H\ Eqs. (1) and (4) of the main 
text, are equivalent in the zero-range limit. 

Columns 2 and 3 of Table U show the coefficients 02,2 
and C3 5 i for the (2, 2) and (3, 1) systems with L u = 1 + 
symmetry and E rel,ni < 21Huj/2. Table [I] shows that ev- 
ery energy shift of the (3, 1) system corresponds to an 
energy shift of the (2, 2) system. However, the (2, 2) sys- 
tem exhibits energy shifts that are not present in the 
(3, 1) system. We find similar results for the n = 4 sys- 
tems with other L u symmetries and for the n = 5 and 
6 systems. As an example, Table [II] summarizes the en- 
ergy shifts for the n = 6 system for the first two non- 
interacting energy manifolds with E rel,ni = 23/kj/2 and 
25hu/2. Columns 3, 4 and 5 of Table HT1 list the 03,3, 04,2 



TABLE I: Perturbative energy shifts for the (2, 2) and (3, 1) 
systems with 1 + symmetry in the weakly- attractive a s < 
regime. The coefficients 02,2 and C3,i are listed for the non- 
interacting energy manifolds ^ rel ' ni — 13ftw/2, 17ftw/2 and 
21hu/2. Within each energy manifold, the coefficients are 
listed in decreasing order. The 03,1 column contains "blank 
rows" to highlight that only a subset of the (3, 1) shifts ap- 
pears in the (2, 2) sequence. The energy shifts for the (2, 2) 
system with E rel > ni < 17huj/2 are taken from Ref. 0. 
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and C5 5 i coefficients for all allowed L u symmetries. It 
can be seen that every energy shift of the (5,1) system 
corresponds to an energy shift of the (4, 2) system. More- 
over, every energy shift of the (4,2) system corresponds 
to an energy shift of the (3, 3) system. However, the (3, 3) 
system exhibits energy shifts that are not present in the 
(4, 2) system, and the (4, 2) system exhibits energy shifts 
that are not present in the (5,1) system. 
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TABLE II: Perturbative energy shifts for the (3, 3), (4, 2) and 
(5, 1) systems in the weakly-attractive a s < regime. The 
coefficients 03,3, C4,2 and cs,i are listed for the two energet- 
ically lowest lying non-interacting energy manifolds. Within 
each energy manifold, the coefficients are listed in decreasing 
order. The C4,2 and cs,i column contain "blank rows" to high- 
light that only a subset of the (3, 3) shifts appears in the (4, 2) 
sequence and that only a subset of the (4, 2) shifts appears in 
the (5, 1) sequence. 
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II. WEAKLY-REPULSIVE a s > REGIME 

In the weakly-repulsive regime, the two-component 
Fermi gas can be thought of as consisting of up to rt2 
composite bosonic molecules and unpaired spin-up and 
spin-down atoms. If no molecules are formed, the per- 
turbative treatment is the same as for the a s < case 
discussed in the previous section. When molecules form, 
unpaired spin- up and spin-down fermions interact, as in 
the a s < case, through the ^-function potential with s- 
wave atom-atom scattering length a s while spin-up atoms 
and molecules (and spin-down atoms and molecules) in- 
teract through a (5-function potential with s-wave atom- 
dimer scattering length a a d [J 0, 0, @|; in addition, the 
atom mass m is replaced by 4m/3 to account for the 
change of the mass of one of the collision partners. If 
two or more molecules form, the molecule-molecule inter- 
action is modeled through the J-function potential with 
5-wave dimer-dimer scattering length a^d and reduced 
mass m [3, 0, 0. Using this effective interaction model, 
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FIG. 1: (Color online) Schematic illustration of effective de- 
scription of two- component Fermi gas in the weakly-repulsive 
small a s /aho regime (a s > 0). Spin- up atoms ("F") are 
depicted by red circles with up-pointing arrow. Spin-down 
atoms ("F'") are depicted by blue circles with down-pointing 
arrow. Composite bosonic molecules or dimers ("D") are de- 
picted by combined blue and red circles. For the (3, 1), (2, 2), 
(4, 1), (3, 2) and (5, 1) systems, the formation of one molecule 
is considered. For the (4, 2) and (3, 3) systems, the formation 
of one molecule and two molecules is considered. 



we have calculated the energy shifts for the n = 4, 5 
and 6 systems within first-order degenerate perturbation 
theory for good portions of the energy spectra. We find 
that the energy shifts of the (n\ + 1, n<i — 1) system form 
a subset of the energy shifts of the (ri\,n<i) system for 
n = 4, 5 and 6. Moreover, the energy shifts of the 
(ni + 2, ri2 — 2) system form a subset of the energy shifts 
of the (ni + 1, Ti2 — 1) system for n = 6. In the following 
we provide a qualitative analysis of the a s > regime 
assuming at least one molecule forms. 

Figure [1] schematically illustrates the effective systems 
encountered for n = 4, 5 and 6 when molecules are 
formed. The (3, 1) system can form one molecule, de- 
noted by "D" (dimer) while the (2, 2) system can form 
one or two molecules. However, as we argue now, only a 
portion of the energy spectrum of the "F+F'+D" sys- 
tem (see Fig. [1]) is degenerate with the energy spec- 
trum of the "2F+D" system. While the wave function 
of the effective (3, 1) system (the "2F+D" system) must 
be anti-symmetric under the exchange of the two spin-up 
fermions, no such constraint exists for the effective (2, 2) 
system (the "F+F'+D" system). This implies that the 
energy spectrum of the effective (3, 1) system containing 
one molecule forms a subset of the energy spectrum of the 
effective (2, 2) system containing one molecule. The anal- 
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ysis of the (4, 1) and (3, 2) systems proceeds similarly. In 
particular, the "3F+D" system contains more symmetry 
constraints than the "2F+F'+D" system. The number 
of atom-dimer interactions are, however, the same. This 
implies that the energy spectrum of the effective (4, 1) 
system containing one molecule forms a subset of the en- 
ergy spectrum of the effective (3,2) system containing 
one molecule. 

For the n = 6 system, the situations where respectively 
one and two bosonic molecules form need to be consid- 
ered separately. We first consider the situation where 
one molecule forms. The degeneracy of a part of the en- 
ergy spectrum of the "3F+F'+D" system with the energy 
spectrum of the "4F+D" system can be understood along 
the same lines as above. For the effective "3F+F'+D" 
and "2F+2F'+D" systems, in contrast, the situation is a 
bit more complicated. Intuitively, one may argue that it 
was shown in Sec. HI that the energy shifts of the "3F+F"' 
system form a subset of the energy shifts of the "2F+2F"' 
system and that the distinguishable "D" particle cannot 
"destroy" the degeneracy. For the (4, 2) and (3, 3) sys- 
tems, two molecules can form. Following the logic above, 
the effective "2F+2D" system contains more constraints 
than the "F+F / +2D" system. This implies that the en- 
ergy spectrum of the effective (4, 2) system containing 
two molecules forms a subset of the energy spectrum of 
the effective (3, 3) system containing two molecules. 

We reiterate that the qualitative pictures discussed 
here have been confirmed by calculating the energy shifts 
for good portions of the energy spectra of the n = 4, 5 
and 6 systems. 



III. UNITARITY 

Our numerical stochastic variational approach for de- 
termining the eigen energies of small trapped few-body 
systems under external spherical symmetric confinement, 
interacting through finite-range Gaussian model poten- 
tials, is discussed in detail in Ref. [7j. For the discussion 
here it is important that this approach yields the ground 
state and excited state energies of trapped few-body sys- 
tems for any s-wave scattering length a s and that the 
accuracy of the eigen energies can be controlled system- 
atically. 

Reference tabulates the extrapolated zero-range en- 
ergies of the (3, 1) and (2, 2) systems described by the 
Hamiltonian H, see Eq. (1) of the main text, at unitar- 
ity. These energies and the associated s%^^ values were 
used to make Fig. 3 of the main text. To exemplar ily 
illustrate the degeneracy of a subset of the (2, 2) energies 
and the (3, 1) energies, Table EI] summarizes the (2,2) 
and (3,1) energies with E rel < 2lhw/2 and L u = 1+ 
symmetry. The second and fourth columns show the cor- 
responding s u ' L,u values, derived from the energies with 
hyperradial quantum number q = 0. It is easily checked 
that the relative energies of the (3, 1) system agree with 
the relative energies of the (2, 2) system to within the 



TABLE III: Extrapolated zero-range energies of the (2, 2) and 
(3, 1) systems with 1 + symmetry at unitarity. The 
values are related to the relative energies at unitarity through 
-^rii,n 2 — (2<? + s ni^ 2 +1)^5 where q = 0, 1, • • • ; the second 
and fourth columns list the s^'/^ values, derived from the 
energies with q = 0. The energies are taken from Ref. |2|. 
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numerical uncertainty of 0.1 % of the extrapolated zero- 
range energies. Similarly good agreement is found for the 
n = 4 systems with other L u symmetries. 

Reference Q also reports the extrapolated zero-range 
energies of the (3, 2) and (4, 1) systems at unitarity; these 
energies have an estimated accuracy of 1 %. For exam- 
ple, the extrapolated zero-range ground state energies of 
the (4, 1) system with L u = 1 + , 2 + and 3 + symmetry de- 
scribed by the Hamiltonian H are 7.155/kj, 7.494fkj, and 
7.502/kj [7[. For comparison, the extrapolated zero-range 
energy of the second excited state of the (3, 2) system 
with L u = 1 + symmetry is 7.158fkj, that of the fourth 
excited state of the (3, 2) system with L u = 2 + symme- 
try is 7.497fk<j, and that of the second excited state of the 
(3, 2) system with L u = 1 + symmetry is 7.504fkj. These 
examples show that the extrapolated zero-range energies 
of the (4, 1) system described by H form, within the nu- 
merical accuracy of Ref. [7[ , a subset of the extrapolated 
zero-range energies of the (3, 2) system described by H. 



IV. APPLICATION OF GROUP THEORETICAL 
FRAMEWORK 

As discussed in the main text, the Hamiltonian H' is 
invariant under the exchange of any pair of particles and 
thus described by the symmetric group S n of degree n. 
In the following, we apply Eqs. (6) and (7) of the main 
text to weakly-interacting atomic Fermi gases (small \a s \) 
with n = 4, 5 and 6 described by H' . 
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We first treat the n = 4 system. The 64 group contains 
five conjugacy classes or irreducible representations that 
are characterized by the Young tableaux [4], [3, 1], [2 2 ] = 
[2, 2], [2, l 2 ] = [2, 1, 1] and [l 4 ] = [1, 1, 1, 1] with degrees 
^4 = 1, d3i — 3, d22 = 2, d 2 n = 3 and dim = 1 0. 
According to Eq. (6) of the main text, the properly anti- 
symmetrized states of the (2, 2) and (3, 1) systems can 
be decomposed as 

[1,1]® [1,1] = [1,1, 1,1] ©[2, 1,1] ©[2, 2] (2) 

and 

[1,1,1]® [1] = [1,1, 1,1] ©[2, 1,1], (3) 

respectively. By inspection, this yields Eq. (7) with n\ = 
ri2 = 2 of the main text, 

[l,l]®[l,l] = ([l,l,l](8)[l])e[2,2]. (4) 

To see how these decompositions relate to the energy 
spectra of the (2, 2) and (3, 1) systems, we determine the 
eigen energies of H f for small \a s \/a\ lo perturbatively. We 
start our discussion by considering the system consist- 
ing of four distinguishable particles with equal masses, 
referred to as the XYZA system. All six two-body in- 
teractions are assumed to be characterized by the same 
s- wave scattering length. Table [TV] lists the coefficients 
c [see Eq. (pQ)], which determine the perturbative energy 
shifts, for each L u symmetry. Column 4 lists the number 
of linearly independent eigen states with zero projection 
quantum number, i.e., the number of eigen states of the 
XYZA system with a given L, II and energy shift (and 
M L = 0). 

To construct eigen functions of four-particle systems 
with a given permutation symmetry, we perform a "post- 
symmetrization" of the linearly independent eigen func- 
tions of the XYZA system. There exist nine physi- 
cal four-body systems with distinct exchange symmetry: 
BBBB, BBBX, BBB'B', BBXY, BBFF, FFXY, FFF'F' 
[referred to earlier as the (2,2) system], FFFX [referred 
to earlier as the (3,1) system], and FFFF. Here, "B" 
stands for "boson" and "F", as before, for "fermion"; 
the particle order does not matter (e.g., the XBBB is 
equivalent to the BBBX system). The number of prop- 
erly symmetrized/anti-symmetrized eigen states for these 
physical systems are listed in columns 5-13 of Table [Vl 
for each L, II and energy shift (Ml = 0). The number 
of linearly independent eigen functions with a given L, II 
and energy shift (and Ml = 0) (see column four of Ta- 
ble HV]) corresponds to the degree of the corresponding 
irreducible representation of the 64 group. 



For the 9ftw/2 energy manifold, e.g., one eigen state 
exists. The "highest symmetry" of this eigen state corre- 
sponds to the BBBB system, which is associated with the 
fully symmetric Young tableau [4] . The BBBB state "au- 
tomatically" also has the proper symmetry of the BBBX, 
BBB'B' and BBXY systems. The three linearly indepen- 
dent states in the llhu/2 energy manifold are associated 
with the Young tableau [3, 1]. Appropriate linear combi- 
nations of these linearly independent functions give the 
properly symmetrized/anti-symmetrized eigen functions 
of the BBBX, BBB'B', BBXY, BBFF and FFXY sys- 
tems. The two linearly independent eigen states in the 
13Huj/2 energy manifold with L u = + and c = 13/2 are 
associated with the Young tableau [2,2]. Appropriate lin- 
ear combinations of these linearly independent functions 
give the properly symmetrized/anti-symmetrized eigen 
functions of the BBB'B', BBXY, FFXY and FFF'F' sys- 
tems. Lastly, the three linearly independent states in the 
13Huj/2 energy manifold with L n = 1 + are associated 
with the Young tableau [2, 1, 1]. Appropriate linear com- 
binations of these linearly independent functions give the 
properly symmetrized/anti-symmetrized eigen functions 
of the BBXY, BBFF, FFXY, FFF'F' and FFFX systems. 

The perturbative calculations summarized above show 
that the eigen states of the FFFX system are, as pre- 
dicted by Eq. (j3]), associated with either the [1,1,1,1] or 
the [2, 1, 1] Young tableaux. In both cases, there exists a 
corresponding eigen state of the FFF'F' system [see also 
Eq. J2J)]. It follows that the eigen energies of the FFFX 
and FFF'F' systems are degenerate. However, the per- 
turbative calculations show — in agreement with Eqs. ([2]) 
and (H]) — that there also exist eigen states of the FFF'F' 
system that are associated with the [2, 2] Young tableau. 
No corresponding eigen states with FFFX symmetry ex- 
ist in this case, explaining the fact that the FFF'F' sys- 
tem supports a set of eigen energies that are absent in 
the FFFX spectrum. In selected cases, we have checked 
that the post-symmetrization approach, applied to the 
unsymmetrized eigenstates of the XYZA system at uni- 
tarity, reproduces the energies of the (2, 2) and (3, 1) sys- 
tems reported in Ref. [7j. Similarly, we tied the energy 
spectra of the (4, 1), (3, 2), (5, 1), (4, 2) and (3, 3) systems 
obtained perturbatively and by the stochastic variational 
approach to the irreducible representations of the sym- 
metric groups 65 and Sq. 



V. ACKNOWLEDGEMENT 

Support by the NSF through grant PHY-0855332 and 
by the ARO is gratefully acknowledged. 



[1] J. von Stecher, C. H. Greene, and D. Blume. Phys. Rev. A 77, 043619 (2008). 

A 76, 053613 (2007). [3] K. M. Daily and D. Blume. Phys. Rev. A 81, 053615 

[2] J. von Stecher, C. H. Greene, and D. Blume. Phys. Rev. (2010). 



5 



TABLE IV: Perturbative analysis of the trapped n = 4 sys- 
tem with various symmetries described by the Hamiltonian 
H' . The first column lists the relative energy of the non- 
interacting energy manifold considered. The perturbative en- 
ergy shifts for a given angular momentum L and parity II are 
reported in column three. The energy shifts are calculated 
without imposing a specific symmetry on the wave functions; 
the corresponding degeneracies are listed in column four (only 
states with projection quantum number Ml — are consid- 
ered). Columns five through 13 list the number of linearly 
independent eigen states with a particular particle symmetry 
that can be constructed from the unsymmetrized states that 
describe the XYZA system. 



£ rel ' ni /(M 


17 


c 


XYZA 


BBBB 


BBBX 


BBB'B' 


BBXY BBFF 


FFXY FFF'F' 


FFFX FFFF 


9/2 


0+ 


12 


1 


1 


1 


1 


i 










11/2 


l - 


8 


3 




1 


1 


9 
z 


i 
i 


i 
i 






13/2 


0+ 


6.5 


2 






1 


1 




i 
i 


i 
i 








9 


3 




1 


1 


o 
Z 


i 
i 


i 
i 










14 


1 


1 


1 


1 


1 
1 












1+ 


4 


3 








1 
1 


i 
i 


o 
z 


i 
i 


1 




2+ 


5 


2 






1 


1 
1 




1 


i 
i 








6 


3 




1 


1 


2 


1 


1 










8 


1 


1 


1 


1 


1 










15/2 


0" 





1 












1 


1 


1 1 




1" 


2.82461 


3 








1 


1 


2 


1 


1 






4.42539 


3 








1 


1 


2 


1 


1 






5 


2 






1 


1 




1 


1 








5.12271 


3 




1 


1 


2 


1 


1 










7.42288 


3 




1 


1 


2 


1 


1 










9.20442 


3 




1 


1 


2 


1 


1 










11 


1 


1 


1 


1 


1 












2~ 


2 


3 








1 


1 


2 


1 


1 






3 


2 






1 


1 




1 


1 








4 


3 




1 


1 


2 


1 


1 








3~ 


3.5 


3 








1 


1 


2 


1 


1 






4 


3 




1 


1 


2 


1 


1 










6 


1 


1 


1 


1 


1 














6.5 


3 




1 


1 


2 


1 


1 







[4] G. V. Skorniakov and K. A. Ter-Martirosian. Sov. Phys. 
JETP 4, 648 (1957) and Zh. Eksp. Teor. Fiz. 31, 775 
(1956). 

[5] D. S. Petrov. Phys. Rev. A 67, 010703(R) (2003). 
[6] D. S. Petrov, C. Salomon, and G. V. Shlyapnikov. Phys. 
Rev. Lett. 93, 090404 (2004). 



[7] D. Rakshit, K. M. Daily, and D. Blume. Phys. Rev. A 85, 
033634 (2012). 

[8] J. J. Sakurai, Modern Quantum Mechanics. Addison Wes- 
ley Longman, Revised Edition (1994). 



